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ABSTRACT 



Context. Complex organic molecules have been detected in massive hot cores for over two decades, and only recently in three hot 
corinos (the inner regions surrounding Sun-like protostars, where the dust temperature exceeds 100 K). Since hot corinos have sizes 
of ~ 100 AU (i.e., of the order of the extent of the Solar System), it is particularly relevant to understand whether they are common and 
to identify the formation route(s) of complex organic molecules. Much has yet to be learned on this topic, since even recent models 
predicted it was not possible to form these molecules in low-mass protostars. 

Aims. We aim to enlarge the number of known hot corinos and carry out a first comparative study with hot cores. The ultimate goal is 
to understand whether complex organic molecules form in the gas phase or on grain surfaces, and what the possible key parameters 
are. 

Methods. We observed millimeter rotational transitions of HCOOH, HCOOCH3 CH3OCH3 CH3CN and C2H5CN in a sample of 
low-mass protostars with the IRAM-30m. Using the rotational diagram method coupled with the information about the sources' 
structure, we calculate the abundances of the observed molecules. To interpret these abundances, we review the proposed formation 
processes of the above molecules. 

Results. We report the detection of HCOOCH3 and/or CH3CN towards NGC1333-IRAS4B and NGC1333-IRAS2A. We find that 
abundance ratios of 0-bearing molecules to methanol or formaldehyde in hot corinos are comparable and about unity, and are relatively 
(depending on how the ratios are determined) higher than those in hot cores and in Galactic center clouds. 

Conclusions. So far, complex organic molecules were detected in all the hot corinos where they were searched for, suggesting that it is 
a common phase for low-mass protostars. While some evidence points to grain-surface synthesis (either in the cold or warm-up phase) 
of these molecules (in particular for HCOOH and HCOOCH3), the present data do not allow us to disregard gas-phase formation. 
More observational, laboratory, and theoretical studies are required to improve our understanding of hot corinos. 

Key words. ISM: abundances — ISM: molecules — stars: formation 
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> 1. Introduction 
• I— I 

Many aspects of the formation of solar-type stars have now 
been elucidated and, even though the details are much debated, 
. . . there exists a widely accepted framewor k for it (see the volume 
Protostars and Planets V, for example. ICeccarelli et al.l 120061: 



IWard-Thompson et al.ll2006HWhite et al.l2006l) . Low-mass stars 
form inside molecular clouds from dense and cold condensa- 
tions, called prestellar cores, which evolve into Class and then 
Class I sources. In the latest phases, the newly born star is sur- 
rounded by a proto-planetary disk, which eventually may form 
planets. 

Somewhat less understood, and therefore more frequently 
debated, is the chemical evolution of matter from the molecu- 
lar cloud to the proto-planetary phase, and then to planets (e.g., 
Ceccarelli et al. 2006). It is now acknowledged that, at the den- 
sities and temperatures typical of the centers of prestellar cores 
(> 10^ cm"^ and < 10 K, respectively), heavy-element bearing 
molecules condense out onto dust grains, forming icy mantles. 
Very likely, hydrogenation and oxidation of small molecules 
and atoms (like CO, O, N, etc.) occur on the grain surfaces, 
so that the mantles end up being composed mainly of water, 
interspersed with traces of formaldehyde, methanol, ammonia, 
and possibly even more complex molecules. As material from 
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the surrounding envelope starts accreting onto the central pro- 
tostar, the increased radiation output heats up the surroundings. 
When and where the dust reaches the appropriate sublimation 
temperature of the ices (which depends on the exact composi- 
tion of these ices), the grain mantles evaporate, injecting their 
components into the gas phase. These components can further 
react to form more complex molecules. The regions where the 
ice mantles sublimate (at Tdust > 100 K) and where the emis- 
sion from complex molecules (whether evaporated from the 
grain m antles or formed in the gas) originates , are called hot 
corinos (' Ceccarellill2004t iBottinelli et alj [2004a). The chemical 
composition of these hot corinos reflects both the heritage from 
the prestellar core phase and the reactions taking place in the 
warm gas. The resul t is a gas rich in complex organic molecules 
jCazaux et alj|2003l) . whose diversity is far from being fully ex- 
plored and understood. It is likely that we have discovered only 
the tip of the iceberg. Moreover, the story does not end here. 
Indeed, the fate of these complex organic molecules is almost 
totally unknown. They may condense onto the grain surfaces 
again during the proto-planetary disk phase. Perhaps they are 
incorporated into the planetesimals forming the building blocks 
of planets, or into comets and asteroids, in which case they may 
end up on newly formed planets as accretion proceeds. This pic- 
ture is, at least partially, supported by the fact that some complex 
organic molecules are found both in hot corinos and in comets 



2 



S. Bottinelli et al.: Hot corinos in NGC1333-rRAS4B and IRAS2A 



(e.g., methyl cyanide CH3CN, methyl formate HCO OCH3, and 
formic acid HCOOH: 'Ehre nfreund & Charnlevll2000 y For these 
reasons, the hot corino phase is not only interesting in itself, but 
it is a critical phase in the process of solar-type star formation. 

In short, the molecular complexity in hot corinos is par- 
ticularly relevant to molecular astrophysics and (exo)planetary 
science, but particularly unknown. Our ignorance can be sum- 
marized by this simple, multiple-part question: which complex 
organic molecules are formed, and why, where, and how does 
this take place? At present we are simply unable to answer 
this question for the following reasons, (i) We do not have a 
census of the complex organic molecules. We suspect that many 
complex molecules are present in the millimeter spectra of hot 
corinos, but the firm identification of a large molecule requires 
several lines from this molecule and can be ham pered by the 
multitude of weak lines in the spectra (e.g., Combe s~et al.lll996t 
ICeccarelli et al.ll20 00c). (ii) We do not know why they form. 
The gas-grain models for hot cores (the high-mass analogs of 
hot corinos) predict that there is not enough time for hot corinos 
to form molecules because the gas falls rapidly towards the 
central forming star before any gas -phase reaction can lead 
to complex organic molecules (e.g.. ISchoier et all l2002h . (Hi) 
We debate whether the obser ved complex molecules reside in 
the passive heated envelope ( Ceccarelli et alj 2000allbl) . in the 
hidden circumstellar disk (j0r gensen et alJ 12005b'), or in the 
interface between the outflow and the envelope ( Chandl er et alj 
l2005h . (iv) Finally, we do not know how complex molecules 
form: in the gas-phase or on the grain surfaces, or both. 
JVIuch of our ignorance stems from the too few hot corinos 
so far studied: IRAS 16293-2422 (hereafter IRAS 16293, 
Cazaux e t al. 2003; Kuanetal. 2004; BottinelH et al. 2004S 
Remijan & Hollis 2006), NGC1333-IRAS4A (hereafter 



IRAS4A, Bottinelli et al. 2004a), NGC1333-IRAS4B (hereafter 
IRAS4B, Sakai et al. 2006), and NGC1333-IRAS2A (hereafter 
IRAS2A. iJ0rgensen et al.ll2005al) . 

In this paper, we present detections of complex organic 
molecules in two Class hot corinos: HCOOCH3 and CH3CN 
in IRAS4B, and CH3CN in IRAS2A. The article is organized 
as follows: details on the choice of sources are given in Sect. 
|2l the observations are described in Sect. [3] and results are pre- 
sented and discussed in Sects. |4] and |5] respectively. Section 
|5] starts with some remarks on the formation routes of com- 
plex organic molecules (Sect. 15.11 based on a summary of 
the possible formation mechanisms of the complex molecules 
HCOOH, HCOOCH3, CH3OCH3, CH3CN, and C2H5CN given 
in Appendix |A|, followed by the analysis of our data (Sects. 
15.21 and 15.31 ), and ends with a comparison of hot corinos with 
hot cores (Sect. lSTil i and with Galactic center clouds (Sect. [53] ). 
Finally, a summary and concluding remarks are given in Sect.|6l 



companion IRAS2B (iLooney et alJl2000t iReipurth et am2002h . 
Two CO bipolar outflows appear to originate within a few arcsec 
of IRAS2A: a highly collimated jet in the east- west direction 
and a large-scale outflow aligned NNE-SSW dSandell et al.l 
119941: iKnee & Sandellll2000l) . 



IRAS 4B and IRAS2A were selected from the sample 
studied in lMaret et~al] (l2004l2005h because they both are good 
hot corino candidates and because, due to their distance and 
luminosity, they are expected to hav e brighter lines com pared to 
other Class sources (e.g., from the lAndre et alj 20001 s a mple) . 
The candidacy relies upon the claim bv iMaret et al.l ( l2004i 
2005) of the presence, in both IRAS2A and IRAS4B, of a warm 
(> 100 K) inner region where grain mantles sublimate. The 
claim is based on the observed jumps in the abundances of 
formaldehyde and methanol: low abundances in the outer, cold 
envelope where formaldehyde and methanol are still frozen onto 
grain surfaces, and high abundances in the inner, warm envelope 
where the heat from the central object causes the desorption 
of formaldehyde and methanol into the gas phase. Indeed, in 
the two sour ces where hot corin os have been detected so far, 
IRAS 16293 ( Cazaux et al.ll2003l) and IRAS4A jBottinelU et all 
2004a), similar jumps in the forma ldehyde and/or methano l 
abundances have also been claimed ( iMaret et alj|2004t l2005h . 
supporting the choice of our targets. 



However, there is a noticeable difference in th e extent of 
the hot corinos (54 to 266 AU, Mar et et al.ll2004 and in flie 
jump sizes in all four sources. Indeed, according to ' Maret et al.l 
(^2005), methanol abundances show jumps of a factor 100 to 
350 in all sources but IR AS4A, where the ju mp is lower than 
about a factor 15, whereas ^Maret et aTl (12004') found formalde- 
hyde abundance jumps of a factor 100 to 6000. Other authors 
have carried out analyses of methanol and formaldehyde in 
low-ni ass protostars (e.g.. ISchoier et al.l 120021: |j0rgensen et al.l 
2005bh. In contrast with what was found bv lMaret et alj ( |2004 . 



J0rgense n et alj (l2005bh claimed that no jump of formaldehyde 
abundance is required to model the line intensities in any of the 
sources exc e pt IRAS 16293 dScho ier et al. 2002). However, like 
Ma ret et alj (l2005h . lJ0rgensen et al.. (2005 b) and 'S choier et all 
(20^ found a meflianol abundance jump in IRAS4B, IRAS2A, 
and IRAS 16293, but not necessarily in IRAS4A. The hot corino 
abundances of formaldehyde and methanol found by the two 
groups a re summarized in Table [T| In the following we will 
adopt the IMaret et all (l2004l2005h fra mework, but we will also 
discuss the results in the light of the |j0rgensen et alj (l2005bl) 
analysis. 



2. Source selection and background 

We observed two Class protostars, IRAS4B and IRAS2A, 
located in the Perseus complex, specifically in the NGC1333 
cloud, whose distance is estimated to be 220 pc dCernisll 19901) . 
IRAS4B belongs to the multiple system IRAS4 and is located 
~ 30" and ~ 17" from the other two components, IRAS4A 
and IRAS4C (iLoonev et al.ll2000l; iReipurth et ani2002l). It may 
be a multiple stellar system itself ( Loonev et al.l " 2000). Both 
IRAS4A and IRAS4B are associated with molecular outflows 
of a dynami cal age of a few thousand years, see n in CO, CS 
(iBlake et all 119951) . and SiO (iLefloch et al.lll998h . IRAS2A is 
part of a protobinary system and is separated by 30" from its 



In summary, the four hot corinos — IRAS 16293, IRAS4A, 
IRAS4B, and IRAS2A — form a sample of interestingly differ- 
ent sources, both because they span an apparent large range of 
formaldehyde and methanol abundances, and because the very 
existence of the abundance jumps is still debated. Hence, even 
though statistically small, this sample will permit a first as- 
sessment of how the presence and abundance of complex or- 
ganic molecules in hot corinos depend on formaldehyde and 
methanol. Indeed, these two molecules are predicted to be 
among the most important parent ones (e.g.,. Casein et al. 19931; 
iRodgers & Charnievll200lL 12003). but it is not always clea r how 
they relate to the complex molecules (e.g.. lHorn et al.ll2004i) . 
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Source 






(HjCO) 




(CH3OH) 




(Le) 


Maret et al. (2004) 


J0i-sensen et al. (2005b)" 


Maret et al. (2005) 


J0rgensen et al. (2()05b)'" 


IRAS 16293-2422 


27 


ixur' 


6x10-"" 


IxK)-' 


IxK)-' 


NGC1333-IRAS4A 


6 


2x10-'* 


3 x 10-' 


< 1 x lO-** 


< 3.5 x lO-' 


NGC1333-IRAS4B 


6 


3x10-'' 


1 X 10-* 


7x10-' 


9.5x10-* 


NGC1333-IRAS2A 


16 


2x10-' 


8x 10-'° 


3x10-' 


1.5x10-' 



" FrQm |j0rgensen et 

^ Note that except in IRAS16293 ISchoier et alJ 2002h. the formaldehyde abundances modeled bv |j0rgensen et in (20053) do not require any jump (see text for further details). 
^ Methanol abundances averaged over A- and E-types. 
From lSchoiefetal] J2003) . 



Table 2, Observed frequencies and targeted molecules. 



Frequency range 


Molecule 




Frequencies 


rm,s (mK) 


rm,s (mK) 


(GHz) 






(GHz) 


1RAS4B 


1RAS2A 


90.07 - 90.32 


HCOOCH3 


A 


90.156, 90.229 


7 


2 




HCOOCH3 


E 


90.145,90.227 


2,7 


2 




HCOOH 




90.164 


2 


2 


98.50-98.75 


HCOOCH3 


A 


98.611,98.682 


7,2 


3 




HCOOCH3 


E 


98.606, 98.711 


7,2 


3 




C2H5CN 




98.610, 98.701 


7,2 


3 


110.29-110.43 


CH,CN 




110.329- 110.383 


5 




223.17-223.42 


CHiOCHj 




223.200 - 223.204 


13 


34 




C2H5CN 




223.385 


13 


34 


226.50 - 227.00 


HCOOCH3 


A 


226.718, 226.778 


11 


15 




HCOOCH3 


E 


226.713, 226.773 


11 


15 


257.35 -257.55 


CH,CN 




257.403 - 257.527 


22 


28 



Note. - A dash indicates that no data were taken in the corresponding 
frequency range. 



3. Observations 

The observations were carried out in June 2003 with 
the 30-meter telescope of the Institut de RadioAstronomie 
Millimetrique (IRAMfl The positions used for pointing were 
(a, (5)(2000) = (03^29"1Z0, 3ri3'09") for IRAS4B, and {a, 
(5)(2000) = (03''28'"55H, 31°1 4'35") for IRAS2A. Based on the 
observations of IRAS 16293 bv lCazaux et al.1 (l2003h . we targeted 
the following complex molecules: methyl formate, HCOOCH3 
(A and E), formic acid, HCOOH, dimethyl ether, CH3OCH3, 
methyl cyanide, CH3CN, and ethyl cyanide, C2H5CN. Different 
telescope settings were used to include as many transitions as 
possible for each molecule (see Table|2|. All lines were observed 
with a low resolution, 1 MHz filter bank of 4 x 256 channels 
split between different receivers, providing a velocity resolution 
of ~ 3, 2, and 1 km s"' at 3, 2, and 1 mm, respectively. Each 
receiver was simultaneously connected to a unit of the autocor- 
relator, with spectral resolutions of 20, 80, or 320 kHz and band- 
widths between 40 and 240 MHz, equivalent to an (unsmoothed) 
velocity resolution of 0.1-0.4 km s"' at 3, 2, and 1 mm. Typical 
system temperatures were 100-200 K, 180-250 K, and 500- 
1500 K, at 3, 2, and 1 mm, respectively. 

Two observation modes were used: position switching with 
the OFF position at an offset of Aa - -100", A6 - -h300", and 
wobbler switching with a 110" throw in azimuth. Pointing and 
focus were regularly checked using planets or strong quasars, 
providing a pointing accuracy of 3". All intensities reported in 
this paper are expressed in units of main-beam brightness tem- 
perature. At 3, 2, and 1 mm, the angular resolution is 24, 16, and 



10" and the main beam efficiency is 76, 69, and 50%, respec- 
tively. Figure[T]shows the obtained spectra. 



4. Results 

Detected transitio ns have been i dentifi ed using the JPL molec 



ular line catalog dPickett et alJ 19981). the Cologne Da tabase 
for Molecular Spectroscopy (" Miiller etaP 1200 iL l2005h . and 
the National Institute of Standards and Technology (NIST) 
Recommended Rest Frequencies for Observed Interstella r 
Molecular Microwave Transitions dLovas & Dragosei |2004|) . 
They are reported in Table[3] We considered only lines with a 3- 
cr detection and a yLSR=7. 0+0.4 km s"' as good identifications. 
We detected two of the five targeted molecules: 4 transitions for 
HCOOCH3-A in IRAS4B and 5 and 6 transitions for CH3CN 
in IRAS4B and IRAS2A, respectively. We also have a possible 
detection for C2H5OH at 226.661 GHz in IRAS4bB Note that 
no observations were made at 1 10 GHz for IRAS2A (see Table 
IIJ, and that the rms reached at 257 GHz for IRAS4B is too high 
to detect the CH3CN transitions at this frequency, if the line 
ratios are similar to those in IRAS4A. 



To derive rotational temperatures and beam-averaged col- 
umn densities (see Table |4|i, we used the rotational diagram 
method (Fig. |2]i. For this analysis, we assumed that the emis- 
sion from a given molecule is unresolved in the smallest beam 
in which a transition of that molecule was detected (flmin, which 
is always > 10" — see cols. 3 and 7 in Table |4]i. This is jus- 
tified by the int erferometric observ ations of I RAS 1 6293-242 2 
(iBottinelli et al.l l2004b; Kua n et al.ll200 4; Chan dler et alll2005h . 
NGC1333-IRAS2A (j0rHensenetaLr|2OO5a), and NGC1333- 
IRAS4A (Bottinelli et al. 2006), which show that the bulk of the 
emission occurs over a region < 1". We therefore corrected the 
integrated intensities observed at lower frequencies for beam di- 
lution with respect to d^in- Then, the hot corino abundances were 
obtained by scaling the beam-averaged colu mn dens ities to the 
sizes of the hot corinos (taken from Table 6 of iMaret et al. 200^ 
^.'25 for IRAS4B and ^.'43 for IRAS2A) and dividing by the hot 
corino H2 column dens ity. The latter was obtai ned from the den- 
sity profiles derived bvl J0rgensen et all (l2002h . based on single- 
dish observations of the dust continuum and radiative transfer 
calculations (see also Table 6 of Maret et al.,..2004.) . This intro- 
duces an uncertainty in the derived abundances, due to the uncer- 
tain sizes and H2 column densities of the hot corinos, caused by 



' IRAM is an international venture supported by INSU/CNRS 
(France), MPG (Germany), and ION (Spain). 



- This transition is not seen in NGC1333-I RAS4A. In this source , 
there is a tentative detection at 90.118 GHz l iBottinellietalJl2004j) 
that has a lower energy (3.5 vs 28 cm-'), but also lower line strength 
(log(I)=-4.8 vs -3.9). 
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NGC1333-IRAS4B 



NGC1333-IRAS4B 



0.04 



0.02 



0.00 



-0.02 




90224 90226 90228 90230 90232 
Rest frequency (MHz) 




NGC1333-IRAS4B 
0.03 r 



0.02 



0.01 



0.00 



-0.01 - 



-0.02 




90154 90156 90158 
Rest frequency (MHz) 



98610 98612 
Rest frequency (MHz) 



NGC1333-IRAS4B 



NGC1333-IRAS2A 



0.06 



0.04 



0.02 - 



0.00 



-0.02 




110340 110350 110360 110370 110380 110390 110400 
Rest frequency (MHz) 




257360 257400 257440 257480 257520 257560 
Rest frequency (MHz) 



Fig. 1. Observed spectra towards rRAS4B and rRAS2A. The rms are 7 mK (top panels) and 5 and 28 mK (bottom panels, left to 
right); the spectral resolutions are 0.2, 0.3, 0.8, and 1.2 km s"' at 90, 98, 110, and 257 GHz, respectively. Detected transitions are 
(from left to right, top to bottom): HCOOCH3-E (80 s - loi), HCOOCH3-A (80 s 25- 624, 836 - 735), CH3CN (6 - 4, 

/sT = 4, 0), and CH3CN (14 - 13, is: = 5, 0). Unlabeled Unes are unidentified.' 



Table 3. Molecular lines detected. 



Molecule 




Transition line 


Frequency 


E/ 






J 7-mbdV 










(MHz) 


(cm-') 


(mK) 


(km s-') 


(Kkms-1) 


IRAS4B 


HCOOCH3 


-A 


72,5 


- 62,4 


90156.5 


13.7 


26 


0.6 


0.017 ± 0.004 






80,8 


- 7o,7 


90229.7 


13.9 


23 


0.5 


0.013 ± 0.004 






83,6 


-73,5 


98611.1 


18.9 


21 


0.7 


0.016 ± 0.004 


HCOOCH3 


-E 


80,8 


- 7o,7 


90227.8 


14.0 


18 


1.8"^ 


0.035 + 0.005 


CHjCN'^ 




64,0 


~ 54 


110349.7 


92.3 


18 


2.3 


0.045 ± 0.009 






63,0 


~ 53 


110364.6 


57.6 


17 


2.8 


0.049 ± 0.011 






62,0 


~ 52,0 


110375.1 


32.8 


24 


2.5 


0.062 ±0.013 






61,0 




110381.5 


17.9 


33 


3.4 


0.119 ±0.014 






60,0 


- 5o_o 


110383.6 


12.9 


42 


3.8 


0.171 ±0.015 


IRAS2A 


CHjCN'^ 




145,0 


- 135,0 


257403.6 


188.5 


169 


3.8 


0.692 ± 0.091 






144,0 


- 134,0 


257448.9 


143.9 


111 


1.6 


0.190 ± 0.046 






143,0 


- 133,0 


257482.7 


109.1 


113 


3.4 


0.413 ± 0.044 






142,0 


- 132,0 


257507.9 


84.3 


145 


2.7 


0.411 ± 0.055 






14,.,, 


- 13,.,, 


257522.5 


69.4 


198 


3.1 


0.662 ± 0.087 






14o.o 


- 13().o 


257527.4 


64.4 


115 


3.8 


0.470 ± 0.068 



" Energy of the upper level of the transition. 
* Width of the observed line. 

^ All the CH3CN lines are unresolved triplets, except at 1 10349.7 MHz, which is an unresolved doublet. The quoted signal is the integral over each 
triplet or doublet. 

This line of the E form of HCOOCH3 has a width about three times as large as that of the corresponding transition of the A form, which could 
be due to the presence of some unknown transition(s) artificially increasing the width. 
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the extrapolation from single-dish observations. However, abun- 
dance ratios are not affected by this uncertainty, unless diiTerent 
molecules are formed in different regions (see discussion in Sect. 

113. 

We make the following remarks about the results we 
obtain. The rotational temperature of HCOOCH3-A is poorly 
constrained due to the absence of points at higher energies. 
More generally, some of the rotational temperatures are below 
100 K. This could be due to non-LTE and/or line opacity effects. 
We find that the derive d CH3CN abund a nce is c onsistent with 
the value obtained by |j0rgensen et al.l (l2005bl) for an inner 
{T > 90 K) region. Upper limits were calculated for the targeted 
molecules that were not detected: HCOOH and C2H5CN in both 
sources, HCOOCH3-A in IRAS2A, and CH3OCH3 in IRAS4B. 
The upper limits for the O-bearing molecules were derived 
using the rotational temperature derived for HCOOCH3-A in 
IRAS4B, whereas the upper limits for N-bearing molecules 
were calculated assuming a rotational temperature equal to that 
of CH3CN in each source. Finally, the CH 3 OCH3 abundance in 
IRAS2A was taken from ' j0rgensen et al.l (12005 al) : the quoted 
abundance (3 x 10"^ in an inner, T > 90 K region) is consistent 
with the upper limit derived from our data (x < 4.2 x 10"^, 
also calculated using the rota tional temp e rature derived for 
HCOOCH3-A in IRAS4B). J0r gensen et al.l (l2005ah also report 
a tentative detection of HCOOCH3, but do not give an estimate 
of the abundance of this molecule. 



5. Discussion 

To investigate whether the data provide information on the for- 
mation of complex molecules, it is necessary to know what the 
possible formation paths are. We therefore report the formation 
reactions found in the literature for the detected molecules in 
Appendix |A1 and we bring up some key points regarding these 
formation routes in Sect. 15.11 Since low-mass protostars were 
thought to have insufficient luminosity to develop a hot core- 
type region, we will then investigate the potential dependence 
of the complex molecules' abundances on the luminosity (Sect. 
I5.2I 1. In Sect. 15.31 we will analyze the hot corino data in view of 
the information given in AppendixlAland look at the differences 
with their massive counterparts and with Galactic center clouds 
in Sects. |574l and |5. 51 respectively. 

5.1. Notes on the formation routes of complex organic 
moiecuies 

For all the complex molecules considered here, we can see from 
Appendix |A] that grain-surface formation is a possible alter- 
native to the "classical" gas-phase formation. In this classical 
view, an estimate of the formation and destruction timescales for 
CH3OCH3 indicates that destruction of this molecule is likely 
compensated by formation mechanisms (Sect. |A.5t . Regarding 
grain-surfaces processes, it is usually implicitly assumed that 
they are taking place during the cold phase preceding the 
warm-up of the dust by the newly born star However, as the 
protostar heats up, it is improbable that the dust temperature 
suddenly jumps from ~10 to ~100 K. Instead, it is more likely 
that the temperature "slowly" rises across the inner envelope , 
leading to a gradual he ating of the grains (e.g.. lVitietaT]l2004l: 
iGarrod & Herbsj 12006'). The effect of this gradual temperature 
change on grain-surface and gas-phas e chemistry can be quite 
substantial. Indeed, the modeling by iGarrod & HerbstI (l2006l) 
shows that gas-phase and grain-surface chemistries are strongly 



coupled during the warm-up phase: molecules formed on the 
grain can evaporate and affect the gas-phase chemistry, whose 
products can re-accrete and in turn change the grain-surface 
chemistry. 

Note that for some molecules, grain-surface formation is not 
only an alternative, but also apparently the only choice. This is 
the case for HCOOCH3 and even more so for HCOOH since 
it has been detected in the ices of st ar forming regions (e.g.. 



ISchutte et alJI 19991: 'Keane et alj|200ll) . as well as of quiescent 
molecular clouds (Knez et alj |2005h . In particular, the detection 
of icy HCOOH in quiescent clouds would indicate not only that 
HCOOH formation occurs on the grain surface, but also that it 
would pre-date the first phase of the star formation mechanism, 
supporting the theory of formation in the cold, rather than 
warm-up phase. Whether the other complex organic molecules 
could also follow the same pattern cannot be commented upon 
since, to our knowledge, no other complex molecule has been 
detected in ices. This is because, unfortunately, the infrared 
spectra of complex organic molecules are not well known in 
this medium. The difficulty in determining the presence of 
these molecules lies in the fact that they may produce only 
slight shifts and broadenings, and that the peak positions 
are characteristic of functional groups, not molecular species 
(Tielens, priv. comm.). Hence, specific identifications are always 
somewhat ambiguous. Overall, it also stands out from Appendix 
|A]that more experiments on thermal surface chemistry at low 
temperatures are needed to check the feasibility of the outlined 
reactions. Cosmic ray processing of icy grain mantles may also 
be able to produce complex molecules, but this process requires 
more quantitative modeling to compare with observations. 



5.2. Luminosity dependence 

As mentioned previously, before the discovery of a hot corino 
around IRAS 16293, it was believed impossible to have chem- 
ically rich regions driven by gas-phase reactions following 
mantle evaporation around low-mass protostars. The argument 
was that, given the low luminosity of these objects, the regions 
where the ices sublimate would be so small that the gas 
crossing-time would be shorter than the time ne eded to form 
complex organic molecules in the gas phase (e.g., ISchoier et al.l 
[20 02) . This is in fact not the case since several hot corinos 
have now been discovered ([Cazaux et al.l 120031 : iBottinelli et al.l 
l2004at |j0rgensen et al.ll2005al this work). However, the ques- 
tion remains regarding the impact of the luminosity on the 
abundances of complex organic molecules in hot corinos. 

Table |5] lists the measured abundances and upper limits for 
the four hot corinos. To remove the uncertainty on the sizes 
of the hot corinos, we choose to look at abundance ratios, in 
particular with respect to formaldehyde and methanol since 
these molecules have been proposed to be the parent molecules 
for complex oxygen-bearing species, if they are formed in 
the gas-phase. These abundance ratios are plotted on Fig. [3] 
as a function of the bolometric luminosity of the low-mass 
protostars. For information, the abundance ratio for CH3CN 
is also plotted (CaHgCN was not included due to the number 
of upper limits), although formaldehyde and methanol are not 
thought to be the parent molecules of nitrogen-bearing species 
(see Sect. |A.4t . Note that the abundance ratios of CH3OCH3 
with respect to both H2CO and CH3OH in IRAS2A seem to 
be "outliers" compared to the other protostars and to other 
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Fig. 2. Rotational diagrams of the detected molecules, corrected for beam dilution at lower frequencies. The arrows show the upper 
limits for undetected transitions. Lines represent the best fit to the data. Error bars are derived assuming a calibration uncertainty of 
10% on top of the statistical error. The excess of emission of the CH3CN transition at 270 K in IRAS2A and 135 K in IRAS4B is 
probably due to contamination from CH3OH 183,16 - 182,17 and unknown line(s), respectively, and are not included in the fits. 

Table 4. Results from the rotational diagrams and upper limits for IRAS4B and IRAS2A. 









IRAS4B 








IRAS2A 




Molecule 






^ beam 








^ ' beam 






(K) 


(") 


(cm--) 




(K) 


(") 


(cm-2) 




HCOOCH3-A 




10 


(4.7 ± 3.6)xlO'' 


(1.1 ± 0.8)xl0-'' 


[38] 


10 


< 2.9 X 10"* 


< 6.7 X 10-' 


HCOOH 


[38] 


22 


< 1.0 X 10'^ 


< l.Ox 10-'^ 


[38] 


25 


< 7.4 X 10'2 


< 1.2 X 10-' 


CH3OCH3 


[38] 


10 


< 6.8 X 10'^^ 


< 1.2x 10-*^ 


[38] 


10 


< 1.8 X 10'" 


< 4.2 X 10-' 


CH3CN 


39 ±3 


22 


(9.7 ± 1.0)xlO" 


(9.5 ±0.2)xl0-* 


87 ± 17 


10 


(3.7 ± 1.0)xl0'2 


(8.7 ± 2.4)xl0-'' 


C2H5CN 


[39] 


10 


< 4.2 X 10'^^ 


< 7.5 X 10-'' 


[87] 


10 


< 4.3 X 10'3 


< 1.0 X 10-' 



Note — When they could not be derived from a rotational diagram, temperatures have been assumed to be equal to that of HCOOCH3 and 
CH3CN for 0-bearing and N-bearing molecules, respectively. In this case, the values are shown in square brackets. 
" Smallest beam size for which a transition was detected (see text for details). 
' Column density averaged over d^i^. 

Abu ndance in the hot corino assuming an H2 column density of A'(H2) = 8.1 x 10^^ cm-- and a hot corino size of 54 AU, i.e., 0.25" at 220pc 
(from lMaret et al.ll2004h . 

Abu ndance in the hot corino assuming an H2 column density of ^(1^2) = 2.1 x 10^^ cm-- and a hot corino size of 94 AU, i.e., 0.43" at 220pc 
(from lMaret et al.ll2004l) . 



O-bearing molecules. However, recall that the CH3OCH3 
abundance was taken from J0rgensen et al. (2005a) where it 
was derived from only one detected transition. Apart from the 
CH3OCH3 points and taking into account the uncertainties 
pertaining to abundance determination, we can see from Fig. [3] 
that the abundance ratios of complex molecules with respect to 
H2CO or CH3OH do not depend on the luminosity, in the range 
~ 5 - 30 Lq. Since the abundances of H2 CO or CH3OH are 
not th emselves a function of luminosity (see lMaret et al.ll2004 
l2005h . then the absolute abundances of the complex species do 
not depend on the luminosity of the protostar. Whatever the 
formation mechanism, either in the gas phase or on the grain 
surfaces, the efficiency in forming complex organic molecules 
is largely constant in the range of studied luminosities. Since 
the luminosity, together with the density, defines the radius at 
which ices sublimate, this also implies that this efficiency is 
rather constant in the inner 200 AUs or so of the studied sources, 
despite the different involved densities (from 10^ to 10'^ cm-^. 



iMaret et al.ll2004 . 

We also investigated the possible dependence of the 
abundance ratios on the ratio of submillimeter to bolometric 
luminosity, Lsmm/Lboi, since it has been suggested as an indicator 
of evolutionary stage. We find that the abundance ra tios do 
not depend on this parameter either. iMaret et aTl (|2004) found 
an apparent anti-correlation between the inner abundance of 
H2CO and Lsmm/Lboi and proposed that this could be explained 
if LsmmA-boi depends on the initial conditions of the protostars 
rather than their evolutionary stage. Indeed, more atomic 
hydrogen is available in less dense (i.e., with a higher Lsmm/Lboi) 
environments, which leads to the formation of more H2CO and 
CH3OH. If we plot the inner CH3OH abundance as a function of 
Lsmm/Lboi, wc also notlcc an apparent anti-correlation. Since the 
abundance ratios are roughly constant with different Lsmm/Lboh 
then the absolute abundances of complex molecules should also 
be anti-correlated with this parameter. Following a similar line 
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of thought as iMaret et al] ( |2004|) and assuming that complex 
O-bearing molecules form on grain surfaces via H, O, OH, 
and/or CH3 additions, this anti-correlation could be indicative 
of these species being more readily available in less dense 
environments. 



5.3. Dependence on methanol and formaldehyde hot corino 
abundances 

Whether we consider gas-phase or grain-surface formation, 
CH3OH and H2CO appear as key molecules: in the first case, 
they have been suggested as parent molecules, in the second, 
they are known mantle constituents. It is therefore interesting 
to investigate the abundance ratios of complex molecules to 
CH3OH and H2CO as a function of CH3OH and H2CO abun- 
dances themselves. We plot these quantities for hot corinos in 
Figs.|4]and|5] respectively. From these figures, we can see that: 

(i) The complex molecules in hot corinos have comparable 
abundance ratios, apparently independent of the CH3OH and 
H2CO abundances. 

(ii) These abundance ratios are close to unity. 

The implications are: 

- In the case of gas-phase formation from methanol or 
formaldehyde: (i) and (ii) mean that, in all the hot cori- 
nos, the formation of complex molecules uses up a sig- 
nificant fraction of the parent molecules sublimated from 
the mantles. However, gas-phase chemistry does not seem 
able to reproduce this behavior For examp le, t h e col- 
lapsing e nvelope model of Rodgers & Char nlevI (120031: 
see also iRodgers & CharnlevI 120011) predicts CH3OCH3 
to CH3C >H abundance ratios of only 10"^ - 10 
Similarlv. lRodgers & CharnlevI ( 1200 lb predicted HCOOCH3 
to CH3OH abundance ratios < 7 x 10"^, and we now 
know that they used too high a formation rate coefficient 
for HC OOCH3, so that t he actual prediction should be even 
smaller dHorn et alj|2004i) . Therefore, either gas-phase mod- 
els are not adequate, or complex molecules are not formed 
in the gas phase. 

- In the case of grain-surface formation: (i) and (ii) show 
that complex molecules are as important mantle constituent 
as CH3OH and H2CO. Observations of soUd HCOOH and 
CH3OH along quiescent lines of si^ ht bv lKnez et al. I (l2005h 
and in protostars ( Keane et a l. 2001) support this idea since 
the quoted HCOOH and CH3OH abundances yield abun- 
dance ratios of order unity, as we find for the hot corinos. 
No observations of other complex molecules in the ices are 
available, but their presence cannot be excluded considering, 
as mentioned in Sect. 15.11 the difficulty of identifying their 
signature in infrared spectra. 

Note that gas-phase reactions proposed in Sect. lA.ll for 
the formation of HCOOH do not involve CH3OH or H2CO. 
Therefore the conclusions mentioned for gas-phase formation do 
not apply to this molecule, whereas the analysis regarding grain- 
surface formation is still valid. 

Using the methanol abundances derived by |j0rgensen et alj 
(l2005bl) does not change the shape of Fig. [4|since thei r values are 
comparable to the ones derived bv'M aret et alj (12005"). However, 
if we take formaldehyde abundances from J0rgensen et al. 
(I2005bl) . we notice that the abundance ratios are larger than those 
plotted in Fig.|5](as expected since lJ0rgensen et alj2005bi do not 



model any H2CO abundance jump), and that the abundance ra- 
tios are scattered by up to two orders of magnitude. 
Overall, in the Maret et al. (2004, 2005) framework, our data are 
consistent with either gas-phase or grain-surface formation of 
complex molecules in hot corinos, with nonetheless more sup- 
port for the later route. However, other scenarios could be con- 
sidered. For example, complex molecules could form in the ISM, 
deplete onto grain mantles during the accretion phase and desorb 
as the protostar heats up its environment. If this formation in the 
ISM were true, we should be able to observe their low-energy 
transitions in dark clouds or in the co ld envelopes of the proto- 
stars. For example, Remiian & HollisI ("2006) observed a transi- 
tion of HCOOH with E„=3.2 K around 1RAS16293B for which 
the spatial distribution is around (and not peaking at) the posi- 
tion determined from the continuum. We cannot therefore rule 
out this theory, but more observations are needed to support it. 

5.4. Comparison hot corinos - hot cores 

Due to the physical differences between hot corinos and hot 
cores, comparing these two types of objects can potentially 
tell us a lot about the parameters influencing the formation 
and evolution of complex molecules. We therefore searched 
the literature for the relevant abundances in hot cores to cal- 
culate abundance ratios in these objects as well. To our knowl- 
edge, there are only four hot cores (Sgr B2(N), Orion hot core, 
Orion compact ridge, and W3(H20)) for which data on all five 
O-bearing molecules studied here (H2CO, CH3OH, HCOOH, 
HCOOCH3, and CH3OCH3) are available, but we can find an- 
other seven sources (Sgr B2(M), G327. 30-0.60, G34.26-i-0.15, 
NGC6334(1), W51 el/e2, G31.41+0.31, and G10.47 -1-0.03) for 
which CH3OH and at least two of the complex molecules 
molecules (HCOOH, HCOOCH3, and CH3OCH3) have mea- 
sured abundances. Data are t aken fro m Gibb et al. (2000a) and 
references therein, and from jkeda et alj (120011) . 

Figures |4| and |5] plot the abundance ratios of massive hot 
cores, as derived from the literature. The first thing to notice is 
that the abundance ratios in the massive hot cores have been 
derived by beam-ave raged column de n sities, except for SgrB2 
(N) and (M), where iNummelin et alj (|2000|) have carried out 
an analysis of the methanol emission similar to what has been 
done for the hot corinos, namely disentangling the hot core from 
the cold outer envelope contribution in the beam. For a given 
comple x mol ecule, if a s ource was listed in both iGibb et al] 
(l2000al) and llkeda et alj (1200 lb . then two abundance ratios 
were derived from the different authors. In this case, the two 
values are reported in in Fig. |4] and linked with a solid line. 
Disturbingly, they can differ by almost two orders of magnitude, 
and the cause of this discrepancy is unclear In Fig. [H we also 
report the results obtained from the hot core + outer envelope 
analysis of CH3OH emission in SgrB2 (N) and (M), to show the 
uncertainty associated with the different methods of abundance 
ratios' determinations. As can be seen, there is only a factor 4 
difference in the abundance ratios if the hot core or the outer 
envelope abundance of methanol is used. 

In principle, the different methods used to derive the 
abundance ratios in hot cores (beam-averaged column density) 
and corinos (full density and chemical structure analysis) 
may lead to such a different result that the comparison may 
be meaningless. However, while this would certainly be the 
case for the absolute values of the abundances, the abundance 
ratios suffer much less from the different methods, as shown 
by the "small" factor 4 in the SgrB2 (N) and (M) sources 



8 S. Bottinelli et al.: Hot corinos in NGC1333-rRAS4B and IRAS2A 

Table 5. Abundances of parent and daughter molecules in our sample of four low-mass protostars. 



Molecule 


IRAS 16293 


IRAS4A 


IRAS4B 


IRAS2A 


Ref. 


H2CO 


1x10"' 


2x10""* 


3x10"" 


2x10"' 


1,2,3 


CH3OH 


1x10-' 


< 1 x IQ-^ 


7x10"' 


3x10"' 


4 


HCOOCH3-A" 


1.7x10-' 


3.4x10"** 


l.lxlO"« 


< 6.7 X 10"' 


5,2,6 


HCOOH 


6.2x10"** 


4.6x10"' 


< 1.0 X 10"* 


< 1.2 X 10"' 


5,2,6 


CH3OCH3 


2.4x10"' 


< 2.8 X 10"** 


< 1.2 X 10"* 


3.0x10"** 


5,2,6,7 


CH3CN 


1.0x10"** 


1.6x10"' 


9.5x10"** 


8.7x10"' 


5,2,6 


C2H5CN 


1.2x10"** 


< 1.2 X 10"' 


< 7.5 X 10"' 


< 1.0 X 10"' 


5,2,6 



Note — Except for H2CO and CH3OH, for which iMaret et alj ( |2004| . [200^ took into account the effects of opacity to derive the abundances of 
these species, it was not possible to determine the optical thickness of the complex molecules' transitions, so that all the complex molecules' 
abundances should be considered as lower limits. 

" We report here the abundances of the A form only of HCOOCH3, for which we have the highest number of rotational diagrams. Note that the 
abundance of the E form is usually very close to that of the A form (from IRAS 16293 and IRAS4A; Cazaux et al. 2003; Bottinelli et al. 2004^, 
so that the total HCOOCH3 abundance would be twice that o f the A form. 

References - ( l)iCeccarelli et al., (i2000a) . (2) iBottineni et alj ( [2004a) . (3) iMaret et al.U2004l) . (4) lMaret et aLH2005l) . (5) ICazaux et"aLH2003h . (6) 
This work. (7') |j0rgensen et al.l ( l2005ah . 
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Fig. 3. The abundances of the observed species normalized to the H2CO (left) and CH3CM'(r\g abundances, plotted as a function 
of bolometric luminosity. Square, triangles, diamonds, and plus signs represent HCOOH, HC OOCH3, CH3OCH3, and CH3CN, 
respectively. The abundance for CH3OCH3 in IRAS2A was taken from |j0rgensen et al] (l2005a) and is likely underestimated (see 
text). Note that we only have an upper limit on the CH3OH abundance in the hot corino of IRAS4A, therefore, we did not plot the 
point corresponding to CH3OCH3 for which only an upper limit is available in that source. Also, we do not show C2H5CN since 
only upper limits are available in three of the four sources. 



dNummelin et al.l [20001) . In practice, the beam-averaged abun- 
dance ratios are very different from the reality only if the spatial 
distributions of the different molecules in the considered source 
are different, and/or if there is a large contribution from the outer 
envelope with respect to the hot core. While a direct measure 
of the molecular emission extent would require carrying out 
interferometric observations, the available data can give a good 
hint on where the emission originates from, by looking at the 
rotational temperature, r,ot. As a matter of fact, the r,ot of the 
hot cores considered in this study (with the possible exception 
of the H2CO and HCOOH) are implying that the emission is 
indeed dominated by the hot cores rather than the cold envelope 



(which may also be due to an observational bias, that is, if the 
observed transitions are probing warm rather than cold gas). In 
this respect, therefore, we think that the comparison between 
hot corinos and hot cores shown in Fig. |4] is reliable. There 
is, however, more uncertainty associated with the hot cores' 
abundance ratios reported in Fig.|5] since H2CO emission could 
have an important contribution from the cold envelope. 

The abundance ratios with respect to CH3OH and H2CO 
were noticed to be roughly constant for hot corinos. Regarding 
hot cores, the ratios seem to decrease with increasing CH3OH 
or H2CO abundance, but the data are also consistent with 
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Fig. 4. Abundance ratios of complex 0-bearing molecules to methanol, plotted as a function of the methanol abundance. Open 
diamonds represent hot corinos. Filled squares represent abundances ratios of hot cores, derived fro m beam-averaged c olumn density 
analysi s. Squares linked by a solid line represent abundance determinations from different authors (iOibb et al.l2000al and references 
therein. Skeda et al]|2001.) . Stars represent the hot cor es of SgrB2 (N) a nd (M), where an analysis of the methanol emission similar 
to what done in the hot corinos has been carried out dNummelin et al.l l2b0d). The dotted lines connect the SgrB2 (N) and (M) hot 
cores (stars) to the squares corresponding to the cold envelopes of these sources. See the text for further details. Note that panel (b) 
represents the total HCOOCH3 abundances, that is twice the abundances of the A form quoted in Table |5] The point corresponding 
to CH3OCH3 in IRAS4A is not represented due to both CH3OCH3 and CH3OH abundances being upper limits in this source. The 
question mark refers to the ratio in IRAS2A and indicates that the CH3OCH3 abundance is likely underestimated in this source (see 
Sectl52]l. 
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Fig. 5. Abundance ratios of complex O-bearing molecules to formaldehyde, plotted as a function of the formaldehyde abundance. 
Open diamonds represent hot corinos. Filled squares represent abundances ratios of hot cores, derived by beam-averaged column 
density analysis. 



a constant ratio with a larger scatter. In any case, it is clear 
from these figures that the abundance ratios with respect to 
methanol in hot cores are lower than in hot corinos by 1-2 
orders of magnitude, whereas abundance ratios with respect 
to formaldehyde are comparable in hot corinos and hot cores. 
We also notice that the HCOOH abundance ratios in hot cores 
are about one order of magnitude lower than the HCOOCH3 
and CH3OCH3 abundance ratios, whereas they are lower by 
only about a factor four in hot corinos. For completeness, we 
considered how Figs. |4] and |5] would change if we were to 
use beam-averaged abundances ratios in hot corinos, as in hot 
cores. In this case, abundance ratios in hot corinos would be 
smaller by a factor < 10: ratios for HCOOCH3 and CH3OCH3 



would become comparable in the two types of objects, while 
HCOOH ratios would still be larger in hot corinos than in hot 
cores. This is very likely due to the relative larger contribution 
of methanol and formaldehyde emission in the cold envelope 
in low mass with respect to high mass protostars. Indeed, high 
energy transitions are more easily detect e d in high than low 
mass protostars (e.g., 'Com ito et alj 120051 : ISchilke et al.l 1200 iL 
1997; Blake et al. 1995), which support the above interpretation. 

Whether considering beam-averaged ratios in all objects or 
not, Fig.|4](and to a lesser extent Fig.jS) show that hot corinos 
are not just scaled versions of hot cores and that in fact, complex 
molecules are relatively more abundant in hot corinos than in 
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hot cores. This conclusion would stil l be valid if we assume 
CH3OH and H2CO abundances from |j0rgensen et al.1 (l2005bh 
since, as mentioned in the previous section, the abundance ratios 
with respect to CH3OH and H2CO would be similar and higher, 
respectively, so in any case, they would be higher than the 
abundance ratios in massive hot cores. A possible explanation 
for the difference between hot cores and hot corinos could lie 
in the grain mantle composition. Indeed, we already know that 
the levels of deuteration differ in these two types of objects, 
more specifically that extreme deuteration occurs in low-m ass 
but not in high-mass protostar s (e.g., Ceccarelli et al] |2006i and 
references therein). Moreover, Boogert et al. (2004) summarize 
the abundances of some mantle constituents (H2O, CO, CO2, 
CH3OH, OCN") in two low-mass and high-mass embedded 
protostars, and large differences can also be seen there. It would 
therefore not be surprising that this would also be the case for 
complex molecules if they formed on the grain surfaces. The 
reason could be, for example, that hot corinos are preceded by 
a longer cold phase during which grain-surface reactions are at 
play, so that complex molecules in hot cores would not have 
time to become an important grain-mantle constituent. 
Alternatively, the difference in abundance ratios could be 
due to a difference in the chemistry, since this depends on 
environmental parameters such as the de nsity and temperature 
of the gas. For example, the model of iRodgers & CharnlevI 
( l200lh shows that the predicted abundances are different at 
100 or 300 K. Also, considering for example the methyl 
fo rmate, the gas-grain mo del in the warm-up phase presented 
in iGarrod & HerbstI (|2006|) shows that the gas-phase formation 
of this molecule is more efficient at low temperatures, and that 
the grain-surface pathway leads to higher abundances the longer 
the dust temperature remains in the 40 - 60 K range. Either 
way, formation of complex molecules during the warm-up 
stage of the protostar would be more efficient in low- than in 
high-mass environments. However, observations performed 
towards Galactic center clouds (see Sect. 15.5b are acting against 
this theory since complex molecules are observed in these 
clouds although they are not subject to a warm-up phase. 
Finally, recall that all the proposed grain-surface formation 
routes of complex molecules involved UV or cosmic -ray 
processing. Hot cores/corinos are shielded from external 
UV radiation fields, but low-mass protostars, unlike mas- 
sive ones, ar e known t o be powerful X-ray sources (e.g., 
iFeigelson &"M ontmerlj Il999l) . X-rays have already been 
proposed to be the re ason for the presence of calcite in the 
low-m ass protostars dChiavassa et al.l 120051: ICeccarelli et al.l 
I2002bl) . In our case. X-rays emitted by Class objects could 
provide the necessary energy to produce a large number of 
radicals and hence a large amount of complex molecules on the 
grain surfaces, thereby explaining the larger abundance ratios 
observed in hot corinos compared to hot cores. 



5.5. Comparison with Gaiactic center clouds 

Hot cores/corinos are not the only objects where complex 
molecules have been observed. Indeed, Requena-Torres et al. 
(I2OO6I) have carried out a survey of complex O-bearing 
molecules in Galactic center (GC) clouds. These clouds 
are known to possess a warm (> 100 K) and not 
too den se (~ 10^ cm~ -^ ) gas (iHiittemeister et al.' 1993'; 
[Rodrig uez-Fernandez et al.l 120001: iCeccarelli et al. 2002a). 
This warm component is probably caused by shoc ks (e.g., 
iFlower et"aLl 119951: iRodrfguez-Fernandez et al.l l2004t) arising 



from cloud-cloud collisions dHiittemeister et all I1993L Il998l) . 
These shocks are thought to be at the origin of the enhanced 
NH3 (Flower et al. 19 93, SiO fe artm-Pintado et al. 1997), and 
C2H5OH dMartm-Pintado et al.ll200ll) abundances. 
Surprisingly enough, considering the very different phys- 
ical environments between GC clouds and hot corinos, 
iRequena-Torres et alj d2006l) found, as we do in the present 
work for hot corinos, that the ratios of the abundances with 
respect to CH3OH (which are comparable to those reported 
in Fig. |4] for hot cores) are approximately constant and do not 
depend on the CH3OH abundance. The authors conclude that the 
gas-phase abundances of the organic molecules they observed 
in GC clouds are likely due to the formation of these molecules 
on the grain surfaces and their release in the gas-phase from 
sputtering/erosion of the grain mantles by the shocks, as it is the 
case for ammonia, silicon oxide, and ethanol. 

Note tha t there is evidence for the pr esence of X-rays in these 
GC clouds dMartm-Pintado et al.ll2000l) . These authors suggest 
that X-rays could contribute to the formation of molecules on 
grain surfaces and evaporate small dust grains. This theory adds 
some support to the possible role played by X-rays in the forma- 
tion of complex organic molecules. 

6. Conclusion 

In this paper, we presented the detections of methyl formate 
and/or methyl cyanide in the low-mass protostars IRAS4B and 
IRAS2A, confirming the presence of a hot corino in their in- 
ner envelope. The conclusions arising from the analysis of these 
observations combined with data on the two other hot corinos 
(IRAS 16293 and IRAS4A) are: 

- Hot corinos are a common phase in the formation of solar- 
type protostars and complex organic molecules are ubiqui- 
tous in Class protostars. 

- The absolute abundances (i.e., with respect to H2) of com- 
plex molecules in hot corinos do not depend on the bolomet- 
ric luminosity, available in less dense environments. 

- Abundance ratios of complex molecules' abundances 
to CH3OH or H2CO abundances (x/x(CH30H) and 
x/ji:(H2C0) are of order unity and do not depend on CH3OH 
or H2CO abundances, indicating that complex molecules 
form on grain surfaces or that gas-phase models have to be 
revised. 

Furthermore, we compared hot corinos with massive hot 
cores. Keeping in mind the different methods of abundance de- 
termination used in the two types of objects, we found that 
x/x(CH30H) and .5t/x(H2CO) in hot cores are relatively lower 
than in hot corinos, showing that complex molecules are rela- 
tively more important in hot corinos, a difference that can be 
explained in the case of grain-surface synthesis. 

Overall, although there is no absolute proof, there is circum- 
stantial evidence for the formation of complex organic molecules 
on grain surfaces, either in the cold phase preceding the begin of 
the star formation process, or in the warm-up phase following 
the birth of the protostar In this scenario, there is also a pos- 
sibility that X-rays emitted by low-mass protostars participate 
in the formation of these complex molecules. Finally, it clearly 
stands out from this work that not only more data (single-dish, 
interferometric, and high-energy transitions) are needed, but also 
that laboratory studies of grain-surface reactions are necessary to 
answer of the question of which, why, where, and how complex 
molecules are formed. 
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Appendix A: Theories of complex molecules 
formation 

In this appendix, we review the different reactions that have 
been proposed for the formation of complex species. These re- 
actions can occur in the gas-phase or on grain surfaces. Gas- 
phase reactions usually involve H2CO or CH3OH as precursors, 
or parents. These are known to form on gr ain surfaces (e.g., 
iTielens & Hagenlll982l: Icharnlev et all 1 19971) and evaporate in 
the gas-phase where they undergo further reaction to form more 
complex, or daughter, molecules. Grain-surface reactions gener- 
ally consist of H or O additions and radical-radical reactions. 

A.1. Formic acid 

Both gas-phase and grain-surface processes have been pro- 
posed for the form ation of HCOOH. In the gas-phase model of 
iLeung et alj (11984 . the precursor ion of HCOOH, HCOOH+, is 
produced via the radiative association 

H2O + HCO+ ^ HCOOH+ + hv (A. 1) 

followed by disso ciative electron re combination to produce 
HCOOH. However. Ilrvine et al.l (1 1990b indicate that they believe 
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that this reaction does not occur due to a competitive exother- 
mic reaction to form H30^ and CO. Instead, they attribute the 
gas-phase formation of HCOOH to occur via an ion-molecule 
reaction followed again by a dissociative electron recombination 

CH4 + ^ HCOOH+ + H (A.2) 
HCOOH+ +e' ^ HCOOH + H. (A.3) 

Reaction iA.2\ has been measured in t he laboratory an d 
found to be quite rapid at low temperature jRowe et al. 
Moreover, the cal culated HCOOH a bundance agrees with the 
values observed bv llrvine et al. l (fT990h in the dark cloud LI 34 N. 

On grain surfaces. iTielens & HapenI (119 82') proposed the for- 
mation of HCOOH through successive additions of H, O, and H 
to solid-state CO: 

CO ±S HCO ±S HCOO ±ii HCOOH (A.4) 

However, radiolysis experiments (si mulating the proces s ing of 
interstellar ices by cosmic rays) by [Hudson & Moord d 19991) 
showed that HCOOH could form in H2O-CO mixed ices via 
the following sequence: 

H + CO ^ HCO (A.5) 
HCO -H OH ^ HCOOH. (A.6) 

The models of lHasegawa & HerbstI (1 1993b are able to reproduce 
the abundances observed in massive hot cores for ages larger 
than 10^-10*' yr Finally, the grain-surface formation of HCOOH 
is suppo rted by the inter ferometric observations of Sgr B2 and 
W51 bv ,LiuetalJ(l2001h . and would also be consistent with ob- 
servations of this molecule in the ices surrounding the massive 
protostar W33A (ISchutte et al.lll999HGibb et alJl2000bl) . 

A.2. Methyl formate 

The commonly accepted formation path for HCOOCH3 starts 
with the reaction between protonated methanol and formalde- 
hyde to form protonated methyl formate and molecular hydro- 
gen: 

[CHaOHzl^ + H2CO ^ [HC(OH)OCH3]+ + H2 (A.7) 

followed by dissociative recombi nation of [HC(O H)OCH3l^ 
with electrons to form HCOOCH3 (Blake et al."l987^. However, 
laboratory and theoretical work by Hornet al. (2004) indicates 
the existence of a very large activation energy for reaction 
( IA.7b . so that the later cannot lead to the formation of proto- 
nated methyl formate. Therefore, the formation of methyl for- 
mate in hot cores cannot occur via this reaction. iHorn et aU 
(j2004|) searched for more favorable transitions between the re- 
actants and products of reaction (IA.7b . but were unsuccessful. 
These authors also investigated reactions involving other abun- 
dant species in hot cores, such as protonated formaldehyde and 
CO. They show that none of the studied processes produces 
enough methyl formate to explain the observed abundances. 
However, they also state that one possibility for producing more 
methyl formate is that formic acid would be synthesized on grain 
surfaces and desorbed into the gas phase, in which case the re- 
action: 

CH+ + HCOOH ^ HC(OH)OCH+ + hv (A. 8) 

would play a signifi cant role. 

The downfall is that lHorn et al] (|2004 find that even if HCOOH 



were injected with an abundance one order of magnitude higher 
than observed in OMC-1, their model still predicts a HCOOCH3 
abundance between one and two orders of magnitude below 
the o bserved value in th is source. Overall, in the light of their 
work, iHorn et aL I (l2004l) conclude that no gas-phase route seem 
able to reproduce the observed abundances of HCOOCH3, 
and hence that this molecule should be produced, at least in 
part, on grain surfaces. As pointed out in Sect. lA.ll there is 
additional evidence in favor of HCOOH being synthesized on 
grain surfaces. 

Two schemes for grain-surface formation of HCOOCH3 
have been proposed, but none of them has undergone laboratory 
investigation yet: 

- Formation fro m precursors CO, O, C, and H landing on grain 



(lHerbstll2005h : 

CO H- H ^ HCO (A.9) 

C ±S CH ±ii CH2 ±ii CH3 ±Q CH3O (A.IO) 

CH3O + HCO ^ HCOOCH3. (A.ll) 



ICharnlev & RodgersI (l2005h mention that many radicals (like 
CH3O and HCO) could form in close proximity via the hot 
secondary electron generated by the passage of a cosmic ray 
through the ice. In this case, it would ensue that radical- 
radical reactio ns such as ( lA.llI ) could occur efficiently. 
- ISorrelll ( 1200 ll) proposed a model in which the photoprocess- 
ing of grain mantles by UV starlight creates a high concen- 
tration of radicals in the bulk interior of mantles. Grain-grain 
collisions then provide excess heat causing radical-radical 
reactions to occur and form large organic molecules. In this 
scheme, HCOOCH3 would be produced from the reaction 
between the carboxyl acid (COOH) and the methyl group 
(CH3) in the following way: 

CO -H OH ^ COOH or HCO + HCO ^ COOH -KSH2) 

CH ±ii CH2 ±ii CH3 (A. 13) 

COOH + CH3 ^ HCOOCH3 . (A. 14) 

However, any grain-surface chemistry preceding the hot 
core/corino phase would occur in a very dense and highly 
visually extinct environment, hence well shielded from UV 
starlight. Therefore, such a UV phot olysis of grain s is un - 
likely to happen, as pointed out by iPeeters et al.l (l2006l) . 
Nevertheless, the radical-radical reaction (IA.14I ) could still 
be a possible formation path for HCOOCH3 via cosmic ray 
processing. 

A.3. Dimethyl ether 

CH3OCH3 was proposed by i Blake et alj ( Il987h to form in the 
gas-phase by methyl cation transfer to methanol, followed by 
electron dissociative recombination: 

CH3OH+ + CH3OH ^ CH3OCH+ + H2O (A. 15) 

CH30CH^ + ^ CH3OCH3 + H. (A. 16) 

iPeeters et al.l (|2006|) claim that their models support such a 
gas-phase route if the methanol abundance is of the order 10"^ 
or more, as it is the case for the OMC-1 hot core (recall that for 
hot co rinos, X(CH30H) < 3 x 10"^). However, Ceccarelli et al] 
(I2OO6I) note that experimental measurements of the neutral 
products of dissociative recombination reactions show that 
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two-body products such as in reaction ( IA.16I ) are often minor 
channels. Therefore, gas-phase formation of CH3OCH3 is 
plausible but not demonstrated. 

On the grains, CH3OCH3 could be produced by a similar 
scheme as for HCOOCH3, that is cosmic ray processing 
followed by the radi cal-radica l reac tion CH3 + CH3O — > 
CH3OCH3 dAUen & Ro binson Il977h . As in the case of 
HCOOCH3, this reaction has not been validated by laboratory 
studies. 



A.4. Methyl and ethyl cyanide 

For completeness, we mention here the possible formation 
routes of these two molecules. However, we will not discuss 
them any further in this study, due to the lack of data on poten- 
tial parents such as NH3. Two substantially different formation 
routes for CH3CN have been proposed in the literature: either in 
the gas-phase or on grain surfaces. 



In the iRodgers & Charnlevl (1200 ll) chemical model of mas- 
sive hot cores, CH3CN is synthesized from NH3 in the fol- 
lowing way: HCN is synthesized from the reaction between 
NH3 and C+ (yielding HCNH+, reaction jAlTT i), followed 
by electron recombination or proton transfer to ammonia 
( IA.I8I 1. CH3CN is then formed from the radiative association 
between the methyl ion and HCN ( IA.19I ), again followed by 
electron recombination ( IA.20b : 

NH3 -H C+ ^ HCNH+ (A. 17) 

HCNH+ -1- e^ HCN -1- H or HCNH+ ^ HCN + Pij8) 
CHJ -I- HCN ^ CH3CNH+ + hv (A. 19) 

CH3CNH+ -H e^ CH3CN + H. (A.20) 



Comparisons between the results from this model and 
observations suggest that a grain-surface formation of 
CH3CN is not required. 

- On grain surfaces, CH3CN can be formed by successive hy- 
drogenation of C2N or by recombination between CN and 
CH3: 

• C2N ±5 HCCN ± S CH2CN ±S CH3C N. This set of 
reactions, used in ICaselU et alJ (Il993h . underestimates 
the CH3CN abundance by a factor about 50 compared to 
the abundance observed in the Orion hot core. This could 
be explained by the fact that HCCN formation is in com- 
petition with C3N formation, and the former is indicated 
by the authors as a less im portant pathway. 

• CN + CH3 ^ CH3CN ("Hasega wa & Herbsll Il993h . 
which yields a CH3CN abundance in good agreement 
with the Orion hot core value for an age of 10^ yr 

Regarding C2H5CN, some studies point towards grain- 
surface formation of this molecule by hydrogenati on of HC3N 
(Blake etal. 1987; Charnley et al. 1992; Caselli etail Il993l) . 
Observations by ILiu & Snvdeii (11999) are consistent with this 
theory. How ever, the abundanc e ratio of CH3CN to C2H5CN 
predicted bv lCaselli et alj d 19931) is at least two orders of magni- 
tude smaller than the ratio observed in hot cores and hot corinos. 
Therefore, it seems unlikely that both C H3CN and C2H5CN 
form via the surface reactions proposed bv lCaselli et al.l (Il993h . 



A.5. Timescales for gas-phase processes 

For the molecules for which no strong argument has been pro- 
posed against gas-phase formation, e.g., CH3OCH3, we can es- 
timate a formation timescale from the slowest rate coefficient of 
the reactions involved. In general, for a reaction between reac- 
tants ri and r2 yielding a product p, we have: 



drip 
IT 



(A.21) 



where kf is the reaction rate coefficient in cm"' s"^. In the event 
that rtrj » rir^ (which is the case for CH3OCH3, as we will 
see further down), we also have n^i + «p ~ «totJ where njot is 
defined as «r, + + rip. Then, 



drip 
dt 



= kf (Mtot - np) rir, 

= -kf rir^ Hp + kf «tot "r. 



rip + kf ntot "r,, 



(A.22) 
(A.23) 
(A.24) 



which has solution of the form np(f) — C\ exp[-f/T/] + C2, 
where Ci and C2 are constants. The formation timescale is 



Tf = (kfti,;) ' 
= (kf Xr^y 



(A.25) 
(A.26) 
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where is the H2 density in cm"^ and x^, is the abundance of 
r2, the least abundant reactant. 

For CH3OCH3, the slowest reaction is ( lA.lSb . which has a 
reaction rate of 1 x lO"'*' cm-' s"' (from the gas-phase chem- 
ical network osu.200^, and (predicted) x^, = -^chsDhJ ~ 
(1 - 5) X 10 '° (from lRodgers & Charnlevl2003h . which is much 
smaller than X,., - xchjOCh, > 1 x 10"^ (from Table |5]l, hence 
justifying the derivation and use of Eq. IA.27I For our sources, 
riH, ~ lO** - 10'' cm-"\ which yields T/ = 6 x 10^ - 3 x lO" yr. 
Regarding the destruction of the studied species in the gas-phase, 
we expect ions such as H30^, HCO^, and H^ to be the main 
destroyers of complex organic molecules during the hot corino 
phase, with H30^ being the most important of the three due 
to the evaporation of water contained in the icy grain mantles 
(ICeccarelli et al.lll996l;lRodgers & Charnlevll2003h . In this case, 
we have: 



drip 

- -kd rip nH30+ 



Td 

ripit) = C3 exp(-f/Trf), 



(A.28) 

(A.29) 
(A.30) 



where kd is the destruction rate, is the destruction timescale, 
and C3 is a constant. Then, 



Td = (^rfnH,o+) 

= {kd "H, -*^H,0+)" 
x-l 



(A.31) 
(A.32) 



10 



10 



Destruction of CH3OCH3 by H3O+ has kd = lO '* cm^^ s ' 
(osu.2005) and xh,o+ = 7 x lO"" - 5 x lO-** 



http://www.physics.ohio-state.edu/~eric/research.html 
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(from Ceccarelli et alj 19961 for the lowest value, and 



iRodgers & CharnlevI l2003l and IStauber et al] l2005l for the 
highest value), which yields Td - 6 - 4 x 10^ yr. This is 
comparable to the mid-range of the formation scale. Moreover, 
as the protostar ages, its luminosity increases, causing the 
expansion of the evaporation front and the replenishment of the 
complex molecules, whether they are mantle species or whether 
they form from freshly evaporated CH3OH and/or H2CO. 
Therefore, the destruction by H30^ is likely compensated by 
the evaporation or formation of more complex molecules. 



